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Abstract 

We present a new adiabatic quantum algorithm for searching over 
structured databases. The new algorithm is optimized using a simpli- 
fied complexity analysis. 

1 Introduction 

Many computational search problems are surprisingly difficult to solve. Quan- 
tum computing is a promising candidate to tackle such difficult tasks. An 
important landmark was achieved by Grover: while an oracle search for 
a marked item out of N classical unstructured items requires, on average, 
0{N) steps, Grover found a quantum algorithm that performs the or- 
acle search in y/N steps. This quadratic speedup leads to the question of 
whether similar speedups can be achieved for other search problems. 

Structured searches are natural extensions of the oracle search. They are 
used when the databases possess some structure. Exploiting the structure of 
the database will increase the performance of a search (classical or quantum). 
Cerf et al [2] were the first to study structured quantum search algorithms. 
Their work is based on quantum circuits. Recently Roland and Cerf [2j gave 
a quasi-adiabatic quantum version of the structured search. 

To understand better the algorithm presented in [21 E] we recall that 
Grover's algorithm uses an iterative improvement strategy. It starts with an 
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equal superposition of all possible states representing items of the database, 
and by iterative application of unitary operations (gates) the initial state 
is rotated towards a state that encodes the solution of the search problem. 
However if the database has some structure of its own other strategies can 
be used. For example, one can use the divide- and- conquer strategy. In this 
strategy one divides the problem into subproblems of manageable size, then 
solves the subproblems. The solutions to the subproblems must then be 
patched back or nested together. For this method to produce good global 
solutions, the subproblems must naturally disjoint, and the division made 
must be an appropriate one (optimal). The optimization results in making 
the division such that the errors made in patching does not offset the gains 
obtained in applying more powerful methods to the subproblems. Roughly 
speaking, the algorithm of |2j uses a combination of the two outlined strate- 
gies. The problem is divided first into two oracle searches, then a third 
global search is used to patch back the two previous searches. 

The present work gives an alternative scheme for dividing the original 
problem into subproblems. As will be seen later, our method allows one to 
relax some of the assumptions made in [2j. 

The paper is organized as follows: in section 2 we present the necessary 
tools for our algorithm, whereas the algorithm itself is presented in section 
3. Using an average complexity analysis for our algorithm, we give in section 
4 the running time for a variety of hard search problems. In the appendix 
we present a generalization of our algorithm to higher nesting levels. 

2 Definitions 

Search problems belong to the family of Constraint Satisfiability Problems 
(CSP). A search problem is defined as the problem of finding a satisfying 
assignment for a set (formula) of constraints acting on n (qu)bits. In general, 
each constraint acts on a number of bits less than or equal to n. The search 
will terminate whenever the "program" finds an assignment, or assignments, 
that satisfy simultaneously all the constraints of the CSP formula. For 
example, the oracle search over an unstructured database corresponds to 
the problem of finding an assignment that satisfies all the constraints of 
the predefined formula, whose constraints are all n bit constraints. Even 
though the number of qubits is fixed, by varying the length or the type of 
constraints the CSP can cover different types of problems. 

For a structured database (problem) we can divide the initial search 
problem of n qubits into subproblems with smaller number of variables. 
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One can choose to test first if an assignment to n A variables out of the 
n variables satisfies all the constraints acting uniquely on this subset of 
variables, labeled A. The result of this test is kept in a register. A similar 
test is then applied to the remaining bits, that we label as subset B. And 
finally one must nest both tests to find a global solution for the initial search 
problem. This is what is meant by a structured or nested search. 

More specifically, assume that the initial search problem admits a num- 
ber Mab of solutions. In a quantum search problem, a solution or a satis- 
fying assignment is a state in the Hilbert space of the n qubits. The set of 
solutions we thereafter denote by 

Mab = {\m t )/i = 1,2, • • • , M AB } . (1) 

All elements of this set satisfy the predefined formula of constraints denoted 
generically by {Cab}- For our patching strategy, the constraints {Cab} are 
classified as 

• {C^}: the set of constraints acting only on the subset A. In general, 
there exist different solutions satisfying all the constrains of this set. 
The set of possible solutions is 

M A = {\mf)/i = 1,2, M A } ■ (2) 

• {Cb}'- the set of constraints acting only on the subset B. Again, 
different solutions can be found to satisfy all the constraints of this 
set. The set of possible solutions is 

Mb = {\mf } /i = 1,2, M B } . (3) 

• | Ca/b } = {Cab} — {Ca} — {Cb} is the set of constraints acting 
simultaneously on A and B . 

Some elements of Ma may not satisfy one or more constraints in jc^/s j, 
those elements can not give global solutions, so they are Not Solutions. The 
set of such elements is denoted by M^ s . On the other hand, the set of 
elements of Ma satisfying all the constraints in {c^/^j, and hence giving 
global solutions, will be denoted by M S A . Therefore 

M A = M S A U M^ s . (4) 

1 This procedure is analogous to the resolution method used in Davis-Putnam procedure 
0. Hence {C A /b} can then be called the reduction of {Cab} by A and B. 
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Similarly we have 

M B = M s b UM b ts . 
So, an element |m,) G Mab can be written as 

\rrii) = \m s A ) <g> \m s B ) , 



(5) 



(6) 



where |m^) and \m B ) are elements of ano - respectively. 

It is evident that for a given state \m s A ) there could be different states 
\m B ) such that their product is an element of M.ab- At some stages, the 
authors of [2] and [3j made the assumption that for a given \mj) there is 
only one \m B ) such that their product is a global solution. Although one 
may argue that such an assumption is justified for hard search problems, 
where there is just one global solution, it is unlikely that it holds for a 
generic search problem. This assumption is not needed for our algorithm, 
and hence our algorithm applies to a broader class of search problems. 

3 Nested algorithm 

The present algorithm consists of two stages. In the first stage (stage /) we 
evolve adiabatically the quantum system from a state which is the ground 
state of a Hamiltonian that can be easily obtained to a state which is a 
product of states of M. A and states of M.b- As previously mentioned, this 
product is not yet a global solution. In stage II a global search, similar to 
what is labeled as stage C in [3], is applied to the output of stage / to rotate 
it towards an element of M.ab- The output of stage II is a global solution 
of the problem. 



In this stage we use the procedure defined in [3]. The Hamiltonian is split 
into two parts, each acting on one of the Hilbert spaces Ha, of dimension 
Na = 2 n - A , and TCb, of dimension N B = 2 ns , spanned by subsets A and B 
respectively. 

The most convenient initial state is the equal superposition of all possible 
pure states of the system 



3.1 Stage I 



*(i = 0)) 
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^oa) ® |*ob) 



(7) 
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where N = 2 n = N a Nb- \^o) is the ground state of the Hamiltonian 

H = (I A - \y A)(^0A\)®lB+lA® (IB- |*0fl)(*0B|) • (8) 

The initial Hamiltonian should evolve adiabatically into a final Hamiltonian 
which has the ground state 

' £ l m ^) "T^P E \ m B) = l*M> ® |*/B> • (9) 



A possible final Hamiltonian is 

H I f = (I A -\¥ fA )(* fA \)®I B +lA®(lB-\*fB)(*fB\) • (10) 

We have chosen Hq and fly so that the evolution is decoupled in TL A and 
TLb- The time dependent Hamiltonian evolving Hq to Hj in a time Tj is 
the linear combination 

H(t) = (1 - s(i)) # + = f(s) H + </(*) (11) 

where s(t) is chosen such that s(0) = and s(Tj) = 1. When the adiabaticity 
condition holds, the initial ground state |^o) will evolve slowly in time Tj to 
the ground state of the final Hamiltonian. The meaning of "slow" evolution 
is quantified by the adiabatic theorem jS], which states that the accuracy 
with which the system remains in its instantaneous ground state is equal to 
the sum of the ratio of the transition matrix element from the ground state, 
with energy Eq, to any other state, with energy Ei, over the fourth power 
of the "radiated" energy Loot = Eq — Ei, i.e. 

IwW^ (12) 



"0i 

This gives a lower bound on the evolution time Tj. 

The eigenvalues and eigenfunctions of the Hamiltonian H(t) can be cal- 
culated analytically. We find that the only nonvanishing transition proba- 
bilities from the ground state are to the two lowest energy levels: 

• E A = Ei + Eq , i.e. the state where the subsystem A is in its first 
excited state, and the subsystem B is in its ground state 

• Eb = Eq + Ei , i.e. the state where the subsystem A is in its ground 
state, and B is in its first excited state. 
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The corresponding transition matrix elements are given by 



{Ei 



where 



.dH. 



\E n 



Mf/Nf / N, 



Mi 



1 = 



£ 2 



A, B 



/ (1 _ 2s )2 + 4^i s (l 



(13) 



(14) 



is the gap (energy difference) between the first excited state and the ground 
state of the i-th subsystem. 

By integrating eq. (|12j) over s using eq. (|13|) the lower bound, or the minimal 
time, will be 

el 

Ti = - I ds 



e Jo 



\ 



=A,B 



(15) 



Note that the parallel evolution, in the two Hilbert spaces TL A and 7i B , of 
stage / is the major feature that distinguishes our algorithm from that of [2]. 
Stage I replaces the two-stages evolution (sequential evolution) used in [2J. A 
priori the sequential evolution seems more efficient, as one can use the result 
of the first stage to eliminate some of the no-good assignments in the second 
stage. However the algorithm used in [2J is not a tree-like search procedure, 
i.e. it does not pick an order in which it instantiates the variables. Hence 
the algorithm of jJJ does not eliminate "bad" trees. Therefore replacing the 
sequential by the parallel evolution will not affect dramatically the time cost 
of the algorithm. On the other hand, as mentioned before , the symmetric 
parallel evolution allows one to analyse a broader set of structured search 
algorithms. 



3.2 Stage II 

The output state of stage I can be written as the sum of a state encoding the 
solution to the search problem and a residual part, which is not a solution: 



IM A M B - M AB 
M A M B 



NS\ 



+ 



where the "solution state" is 



1^ 



1 
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(16) 
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and we have used the fact that MJmJ = Mab- 

At this stage one is inclined to start again an adiabatic evolution from 
a Hamiltonian whose ground state is I^/a) <8> e.g. 

fli = l-|*/A>(*/A|®|*/B><*/B| (18) 
to a final Hamiltonian whose ground state is \^ s ), e.g. 

H f = l-\^ s ){^ s \ , (19) 



in time 0(\J MaMb/Mab )■ Unfortunately, this is not possible since the 
initial Hamiltonian Hi is not accessible. Following [3] a global search can 
be applied to nest between the two subsets A and B. The global search is 
achieved through the following procedure 

• Stage / is approximated by an evolution operator U such that ® 
|*/ B > ~ U\^ 0A ) ® |*ob}- 

• Hence Hi ~ UHqJJ^, note that ffj is "replaced" by Ho which is easily 
accessible. 

• The adiabatic evolution from Hi to Hf is implemented on quantum 
circuits. The continuous evolution is replaced by a "step evolution" 
over intervals of time. In each step, at a given time t, we use the 
approximation 

e -iH t t _ jje-iHotjjt ( 2Q ) 

i.e. a backword evolution in time Tj, an application of exp(-iHot), 
then a forward evolution in time Tj. This step requires a time 0(Tj). 

The number of steps is chosen to minimize the error involved in the dis- 
cretization needed during the global evolution. The number of iteration is 
at least (see 0) of the order of y/ MaMb / Mab ■ So that the total running 
time needed to get a global solution, is proportional to 



This expression is not yet the final answer, since there is still the problem 
of determining the values of Ma,b,ab- This will be discussed in the next 
section. 
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4 Complexity analysis 



The total running time T is a function of Ma, Mb, Mab, n and x = ua/ti, 
among which n is the only input parameter, and x will be chosen to minimize 
T. Hence we are left with the three parameters Ma, Mb and Mab- A pri- 
ori, these parameters are problem dependent and are difficult to determine. 
However, since we are interested in getting general results which depend as 
little as possible on the details of the problem, we can approximate these 
quantities using a complexity analysis. 

Computational complexity theory studies the quantitative laws which 
govern computing. It seeks a comprehensive classification of problems by 
their intrinsic difficulty and an understanding of what makes these problems 
hard to compute. The simplified average complexity analysis derived in 
E| can be adapted for our generic algorithm. The major simplification 
made in j^j is the approximation of independent no-good assignments 2 . This 
allows for a relatively easy combinatoric analysis. Following |2j the unknown 
parameters are estimated by 

Mj « 2 n i- na{ - n i/^ k for j = A,B,AB with n AB = n , (22) 

where k is the number of variables per constraint, assuming a constant 
length constraint, a represents the average difficulty of the problem, it 
characterizes the average number of no-good ground instances per variable. 
Hardest problems are found near a critical 3 value a = a c = 1. The critical 
a is obtained from the number of solutions 

M AB ~ N 1 -* 

for a < 1, Mab is large, i.e. the problem is under-constrained and easy to 
solve, while for a > 1, Mab < 1 5 i-e. the problem is over-constrained and 
probably has no solution. Finally for a = 1 we have one solution, this is the 
definition of the critical value. 

Therefore, the optimization procedure reduces to a minimization of the 
running time as a function of the ratio x for different values of a and k. 

2 For most CSP there are classes of no-good assignments that can be deduced from each 
other. For example, one can use proposional reasoning to generate some no-goods from 
others, see for example chapter two of |S]. 

3 This critical behavior is similar to a phase transition in condensed matter physics. 
Similar phase transition features were observed for random K-SAT problems |§lll(J|. 
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4.1 Approximate and exact scalings 

The running time T can not be calculated analytically. However by the 
following simple arguments it is possible to give an estimate for T: 

• We first recognize that our algorithm is inefficient if Mi/Ni ~ 1. This 
can be easily seen by remembering that iVj is the number of possible 
assignments in space A or B and Mi is the number of solution in that 
space. For Mj ~ Ni any assignment will become a solution, then, 
one step (picking up any state) is sufficient to find a solution in the 
subspace i. This is reflected through the limit 




& = ' / _ 1 _ i _> o when M< -► N, 



Therefore, the first stage does not eliminate any no-good assignment 
and consequently the second stage will take VN steps to give a so- 
lution, i.e. as if the search is an oracle search. But what does this 
means? 

Using the average complexity analysis of the previous section we write: 



M A /N A = 2 



-nax 



Taking Ma/Na ~ 1 implies that nax k ~ 1 or k ~ ln(na)/ln(l/x), on 
the other hand we expect that optimization will give a value of x not 
far from 1/2, hence for large values of n and near the critical a = 1 
our algorithm breaks down for k ~ ln(n). But the larger the k is, the 
closer we are to the unstructured search where each constraint is an n 
qubit constraint. 

• From the previous remark we deduce that our algorithm is effective as 
long as Mi/Ni <C 1. 

• The running time of stage /, given in eq. (|15j) . exhibits a near singu- 
lar behavior for Mi/Ni <C 1. This near singularity occurs when the 
integration variable s is close to the value for which f{s) = g(s), i.e. 
s = 1/2. 

These remarks imply that, for a given x, Tj can be approximated using the 
near singular behavior to obtain 



1 If x = or 1, the first stage is again useless as there is no division of the initial space. 
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Figure 1: The numerical value of the running time T as a function of x, 
for (k = 2, a = a c = l,n = 32), compared to the analytical estimate of T. 
Similar results are obtained for different values of k. 



In other words, the time Tj scales with the square root of the largest value 
of the dimension of each Hilbert space divided by the number of possible 
solutions in that space. Multiplying by the number of iterations of stage II, 
the total running time reads 



ti \ L ( N aMb N B M A \ 

Using the average complexity analysis we get 

T(x) ~ y^Max (N<*-*(.i-x) k , N a - axk ^ . (25) 

In figure ^ we compare this approximate formula with the time calculated 
numerically. The two results match extremely well for all values of x. 
The time T is then optimal for jV"-"(i-z) fc = N a-ax k ^ i e x = which 
is compatible with the symmetric nature of the algorithm. This result agrees 
with the numerical results presented in figure |2J where x = 1/2 is found to 
give the optimal time for different values of k. The plot also illustrates that 
the larger k is the closer we are to an oracle type search, which is what we 
predicted before. Hence the optimal time is given by 

Tmin = T(X = 1/2) ~ iV«/2-«/2 fc+1 „ 2 na(l/2-l/2*+i) (2g) 

We conclude that the parameter of interest in these considerations is na. In 
figure 01 we plot the running time as a function of x for different values of 
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Figure 2: The running time T as a function of x, for [k = 2, 3, 5, a = a c = 
l,n = 32). The optimal time is attained at x = 1/2. The horizontal line is 
the y/~N scaling. It is evident that the larger k is the closer we are to the 
horizontal line, i.e. to the oracle search. 

a and n keeping their product constant. The fact that the different plots 
coincide strongly supports the scaling obtained from the analytical estimate. 
Another important aspect is the (sub)-exponential growth of T as a function 
of a. This can be seen either from eq. (|26|) or from the "critical" behavior 
depicted in figure EJ 

Finally we consider the special case (k = 2, a = 1). In this case the 
running time is 0(iV 3 / 8 ), which is better than the classical running time 
O(iV - 5 ) but less efficient than the result obtained using the algorithm of 
Cerf et al. which gives O(iV ' 31 ). The O(iV ' 375 ) scaling is compared with 
the numerical result in figure El 

5 Conclusions 

We have presented a new adiabatic quantum algorithm for searches over 
structured databases. Our algorithm is constructed with minimum assump- 
tions about the nature of the database and the specificity of the search 
problem. Moreover, it results in a significant potential speedup over its 
classical counterpart. 

The new algorithm is in fact "quasi" -adiabatic, it requires the use of 
quantum circuits at some stages. An interesting topic for future work would 
be to construct a "pure" adiabatic quantum algorithm. Another important 
issue is to explore the possibility of using smart classical strategies such 
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Figure 3: The running time T as a function of x, for (k = 3,a = 1.054, n = 
32), (k = 3, a = 1.25, n = 27), such that na is approximately constant. This 
supports the effective parameterization in terms of na. 
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Figure 4: The running time T as a function of x, for {k = 5, a = 
0.9, 1, 1.054, n = 32). The slightest variation of a affects dramatically the 
running time, which reflects the critical behavior predicted by the complex- 
ity analysis. 
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Figure 5: The running time T as a function of N, for {k = 2,x = 0.5). 
The numerical result matches the approximate analytical scaling: T ~ 
O(iV - 375 ). 



as back tracking and constraint propagation, and to see the effect of such 
strategies on quantum interference and coherence. These and other related 
issues will be addressed elsewhere. 



A Multi-nesting verses mult i- part it ion 



The nesting procedure can be applied to subsets A and B, resulting in 
higher level nesting or multi-nesting [2]. In multi- nesting stages / and II 
are applied to subsets A and B separately. Then stage 77 is applied to nest 
A and B. Then the time Tj is replaced by a shorter time. More nesting can 
be used to enhance the effectiveness of the structured search. 

Another alternative is to modify stage / by splitting the n variables into 
more than two subset then 



1 r 1 

Tj = - ds 
e Jo 



\ 



E 
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(27) 



i=A,B,C,D, 

and hence applying the procedure of stage II we get 



Ti 



n M t 

i=A,B,C,D,... 

Mab 



T(M A ,M B ,M c ,...,M AB ,n,n A /n) 



(28) 



An average complexity analysis can then be applied to obtain an estimate 
for T. 
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